Response of liquid xenon to Compton electrons down to 1.5 keV 
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The response of liquid xenon to low-energy electronic recoils is relevant in the search for dark- 
matter candidates which interact predominantly with atomic electrons in the medium, such as 
axions or axion-like particles, as opposed to weakly interacting massive particles which are predicted 
to scatter with atomic nuclei. Recently, liquid-xenon scintillation light has been observed from 
electronic recoils down to 2.1 keV, but without applied electric fields that are used in most xenon 
dark matter searches. Applied electric fields can reduce the scintillation yield by hindering the 
electron-ion recombination process that produces most of the scintillation photons. We present new 
results of liquid xenon's scintillation emission in response to electronic recoils as low as 1.5 keV, 
with and without an applied electric field. At zero field, a reduced scintillation output per unit 
deposited energy is observed below lOkeV, dropping to nearly 40% of its value at higher energies. 
With an applied electric field of 450V/cm, we observe a reduction of the scintillation output to 
about 75% relative to the value at zero field. We see no significant energy dependence of this value 
between 1.5 keV and 7.8 keV. With these results, we estimate the electronic-recoil energy thresholds 
of ZEPLIN-III, XENON10, XENON100, and XMASS to be 2.4keV, 1.8keV, 1.7keV, and 1.1 keV, 
respectively, validating their excellent sensitivity to low-energy electronic recoils. 



PACS numbers: 95.35.+d, 14.80.Va, 29.40.Mc, 78.70.-g, 61.25.Bi 
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I. INTRODUCTION 

Liquid xenon (LXc) provides an ideal detection 
medium for Weakly Interacting Massive Particles 
(WIMPs), which are promising and testable candidates 
for cold dark matter in the Milky Way [1]. As a noble 
liquid, xenon is both a very good scintillator and ionizer 
in response to the passage of radiation [2]. The possi- 
bility to detect charge and light signals after a WIMP 
scatters on a xenon nucleus, along with the relative ease 
to scale-up to large masses and its self-shielding proper- 
ties (e.g. high stopping power for penetrating radiation) 
has made LXe the WIMP target of choice for many at- 
tempts to directly observe a WIMP-induced signal in the 
laboratory [3-7]. 

The expected signature of WIMP interactions in a 
LXe detector are low-energy depositions by the recoil- 
ing xenon nuclei up to several tens of keV. A fraction 
of the energy of the recoiling nucleus is transferred to 
electronic excitations of the medium and is observable 
as ionization and scintillation. The energy dependence 
of the ionization and scintillation signals in response to 
nuclear recoils has been investigated by multiple groups 
below lOkeV [8-11], and most recently down to 3keV 
[12]. In contrast, until recently [13, 14], studies of LXc' s 



'Present address: Institut fur Physik, Johannes Gutenberg Uni- 

versitat Mainz, 55099 Mainz, Germany 

' Corresponding author: aaronmlSphysik .uzh . ch 

* Present address: Max-Plank-Institut fur Kcrnphysik, 

Saupfercheckwcg 1, 69117 Heidelberg, Germany 

§ Present address: Albert Einstein Center for Fundamental Physics, 

University of Bern, Sidlerstr. 5, 3012 Bern, Switzerland 



response to low-energy electronic recoils have not been 
pursued with the same voracity that has accompanied the 
measurements of nuclear recoils. This is not surprising, 
as electronic recoils present only background in searches 
for WIMPs. However, dark matter searches that focus 
on WIMPs can also have sensitivity to non-WIMP dark- 
matter candidates that may interact with electrons. For 
example, axion-like pscudoscalars, which couple to two 
photons, could ionize an atom via the axioclectric effect, 
which is analogous to the photoelectric process [15, 16]. 
This process has been considered as a possible dark mat- 
ter interpretation of the annual modulation signal ob- 
served by the DAM A/LIBRA experiment in the region 
around 2-5 keV [17, 18]. Part of the parameter space 
consistent with this interpretation has been excluded by 
the CoGeNT and CDMS-II experiments [16, 19]. LXc 
could, in principle, have sensitivity to much (if not all) 
of the remaining parameter space, but a measurement 
of LXe's response to low-energy electronic recoils under 
the conditions present in most LXe dark matter searches 
(i.e. applied electric fields) has yet to be achieved. 

We report here a study of LXe's scintillation yield to 
electronic recoils as low as 1.5 keV resulting from low- 
angle Compton scatters of 662 keV 7 rays. For this pur- 
pose, a small LXe cylindrical cell is irradiated with these 
7 rays, a fraction of which are absorbed by a Nal scin- 
tillating crystal placed at a set of chosen angles relative 
to the original direction of the 7 rays, thus kinemati- 
cally determining the energy of the recoiling electrons. 
Recently, another group has used a similar technique to 
measure the LXe scintillation yield down to 2.1 keV [14]. 
We additionally measure the scintillation quenching in- 
duced by the application of a uniform electric held to the 
LXc, which is essential to infer the true energy thresholds 
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FIG. 1: Schematic top-view of the experimental setup. The 662 keV 7 rays are collimated twice: first as they leave the 137 Cs 
source, and second after they scatter in the LXe volume. The Pb channel from LXe to Nal is also covered on top and bottom 
(not shown). The scattering angle, 6, is varied from 4.25° to 34.5°. 



of existing LXe dark matter searches, many of which ap- 
ply fields of similar strength [3-6] . This study represents 
the first observation of scintillation signals (both with 
and without applied fields) in this energy range. The pa- 
per is organized as follows: in Section II we describe the 
experimental methods used for the Compton scattering 
experiment, namely the LXe cell and the Nal detector. 
Section III describes the Monte Carlo methods used in 
simulations of the setup, and in Section IV we present the 
data analysis, including comparison with detailed Monte 
Carlo simulations, and give the results of our measure- 
ments. In Section V we present a summary of our main 
findings, as well as a discussion and implications of the 
results for dark matter searches. 



II. EXPERIMENTAL METHODS 

The Compton-scatter setup consists of a collimated 
137 Cs source, a small LXe scintillation cell, and a Nal 
scintillating crystal, shown schematically in Figure 1. 
The 17.3 MBq 137 Cs source emits 662 keV 7 rays and is 
encased in a lead block with a small cylindrical open- 
ing, 0.6 cm in diameter and 5 cm long, that acts 
collimator. Monte Carlo (MC) simulations of this source 
show that the resulting beam from the collimator has 
a 1 cr angular spread of 1.6°. The LXe cell, which is 
described in detail in [13, 20], consists of a cylinder of 
LXe, 4.5 cm tall and 3.5 cm diameter, viewed on top and 
bottom by two 2" -diameter Hamamatsu R6041 photo- 
multiplier tubes (PMTs), and surrounded by a polyte- 
trafluoroethylene (PTFE) shell. The PTFE acts as an 
efficient light reflector [21] which permits photons hit- 
ting the detector walls to still be detected in the PMTs. 
Three flat grid electrodes, located at 0.5cm (cathode), 
3.5cm (gate), and 4cm (anode) above the bottom pho- 
tocathode, intersect the LXe cylinder and are used to 
apply static electric fields across the volume. In order 
to maximize the efficiency for detecting scintillation pho- 
tons, LXe is filled fully from the bottom PMT to the top 
PMT, producing a single-phase detector. This contrasts 
with most LXe dark matter detectors which use a dual- 
phase design in order to also detect very small ionization 



signals [22]; the scintillation signal in the present detec- 
tor is reduced by ^40% when the liquid-gas interface is 
lowered below the top PMT. The PMT photocathodes 
are held at ground potential, with positive high voltage 
applied to their anodes. Throughout the run, the LXe is 
continuously recirculated and purified through a SAES 
Monotorr hot getter, in order to remove any impurities 
that may enter the liquid. The Nal detector is a Saint- 
Gobain model 3M3/3, which is a fully integrated crystal 
and PMT. The Nal crystal itself is a cylinder, 7.6 cm in 
diameter and in 7.6 cm height. 



The opening of the source collimator is placed initially 
70 cm from the center of the LXe cell. For a subset of the 
scattering angles (4.25°, 5.25°, and 8.5°) this distance is 
reduced to 28 cm (the minimum allowed given the detec- 
tor components) in order to minimize the beam's spot 
size within the LXe volume. A distance of ~1 m is cho- 
sen for the Nal position as a compromise between event 
rate, which decreases with larger separations, and an- 
gular systematics (see Section III), which improves with 
increased separation. The three components are aligned 
using a goniometer with 0.25° tick marks; this tick-mark 
width is taken to be the la accuracy (±0.125°) of the 
geometrical alignment and is included as a systematic 
uncertainty in the analysis (see Section IV). The pre- 
cision with which a scattering angle can be reproduced 
is better than the spacing between adjacent tick marks, 
and therefore associating this width as a ler uncertainty 
is conservative. Unless otherwise specified, reported scat- 
tering angles refer to the angle formed by the collimated 
beam with the centers of the detector components. After 
scattering in the LXe cell, the 7 rays are further col- 
limated on their way to the Nal detector by means of 
a lead channel with a 3 cm circular aperture at its en- 
trance (LXe side), which then widens to encompass the 
Nal crystal and PMT (see Figure 1). Data are collected 
at central scattering angles of 4.25°, 5.25°, 6.25°, 8.5°, 
16.25°, and 34.5°. These correspond to expected elec- 
tron energies of 2.35 keV, 3.57 keV, 5.05 keV, 9.28 keV, 
32.5 keV, and 123 keV, respectively, when applying the 
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well known Compton scatter formula, 
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where E er is the energy of the recoiling electron, E 1 is 
the initial energy of the incident 7 ray, m e is the mass of 
the electron, and 9 is the scattering angle. However, as 
will be shown in Section III, the finite size of the detector 
components lead to peak recoil energies that differ from 
these expectations. 

All three PMT signals — two from the LXe and one 
from the Nal — are read out directly, without amplifica- 
tion. The signals are each split with a CAEN N625 lin- 
ear Fan-out. One triplet of signals from the Fan-out are 
passed to the trigger and timing system, while the other 
triplet is fed directly into an Acqiris DC436 lOOMS/s 
waveform digitizer. The signals going to the trigger sys- 
tem are first passed to a CAEN N840 leading edge dis- 
criminator, where they are converted to logical signals. 
Individual trigger thresholds for the two LXe channels are 
set at ~1 photoelectron (PE). The global trigger condi- 
tion requires time coincidence between the Nal signal and 
either of the two LXe signals. In this way, the system 
utilizes a two-fold PMT coincidence requirement globally, 
but allows the LXe itself to trigger with only a single 
channel, thus optimizing the LXe scintillation detection 
efficiency. The trigger detection efficiency is directly mea- 
sured with a 22 Na source, similar to the technique used in 
[12]. The resulting detection efficiency curve is shown in 
Figure 2, reaching unity by ~2PE. The smallest signals 
considered for analysis (Section IV) are 3PE. 
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FIG. 2: The measured hardware detection efficiency of the 
experimental setup, as a function of the measured scintillation 
size in the LXe. The lowest analysis threshold considered is 
2 PE, where the detection efficiency is already unity. 



The logical signals from the LXe and Nal are addi- 
tionally connected to the "start" and "stop" inputs of 
an Ortec 566 Time to Amplitude Converter (TAC). Its 
output is digitized, along with the three PMT signals, to 
measure the Time of Flight (ToF) of the 7 ray and further 
constrain the time coincidence of the LXe and Nal signals 
offline. The TAC module has an intrinsic timing resolu- 
tion of better than 5 ps, however, the true precision of the 
ToF measurement is limited by the rise time of the Nal 



scintillation pulse, which is of O(ns). The ToF measure- 
ment method is calibrated with the use of a 22 Na source 
placed between the LXe and Nal. This source undergoes 
f3 + decay, eventually emitting two 511 keV photons (in 
opposite directions) when the positron loses energy and 
annihilates with an electron; a signal in both detectors 
physically indicates ToF= and can be shifted with vari- 
able delay generator to calibrate the full range of the TAC 
system. 
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FIG. 3: (color online) Measured spectra in the LXe from the 
decay of 57 Co (122 keV and 136 keV) and 83m Kr (32.1 keV and 
9.4 keV). The feature in the 57 Co spectrum at and below ~500 
PE is the Compton continuum from this source. 



The PMT gains and light-yield stability are calibrated 
weekly, throughout the duration of the data taking. 
The LXe PMT gains are measured by means of a low- 
intensity, pulsed, blue light emitting diode placed in 
the LXe but outside of the PTFE shell. This place- 
ment is chosen in order to help diffuse the light reach- 
ing the PMTs, thereby illuminating all regions of the 
photocathodes. The procedure is similar to that used 
in experiments such as Borexino [23] and XENON100 
[24]. With the measured gains, all PMT signals from the 
LXe are converted to units of PE. The light yield sta- 
bility is monitored using two radioactive isotopes, Co 
(external, 122 keV and 136 keV 7 rays) and 83m Kr (in- 
ternal, 32.1 keV and 9.4 keV transitions, mostly conver- 
sion electrons). The spectra from one set of calibra- 
tions with these sources are seen in Figure 3. The two 
7 rays emitted by the 57 Co are not separately distin- 
guishable, and instead produce a broadened peak whose 
mean value inside the LXe, determined by MC simula- 
tions, is 126.1 keV. The LY obtained from this source is 
(13.05 ± 0.04)PE/keV. The use of 83m Kr in this type 
of particle detector is a relatively new procedure, and is 
described in detail in [13, 25]. 

For most of the measurements, the three grid elec- 
trodes in the LXe are grounded, along with the PMT 
photocathodes, ensuring that there are no electric fields 
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within the LXe volume. However, for three of the central 
scattering angles studied, data are also collected while 
applying high voltage (HV) to these electrodes, in order 
to measure the light yield's field dependence. For this 
purpose, the cathode, gate, and anode grids are biased 
at -265 V, -530 V, and -2120 V, respectively. Due to 
effects of field leakage through the grid wires, the field 
in the main LXe volume (between the cathode and the 
gate, where the 7-ray beam is centered) is less than what 
would be naively determined by treating the electrodes as 
infinite solid plates. In order to study the electric fields 
in all regions of the detector, its geometry and electrode 
voltages are simulated with the comsol Multiphasics® 
software (version 4.3) [26]. These simulations show that 
the volume-averaged field strength in the main LXe re- 
gion is |E| = (450±8)V/cm [27]. 



III. MONTE-CARLO SIMULATIONS 

Given the small scattering angles involved in this 
study, a detailed Monte Carlo (MC) simulation of the 
setup is essential to understand its systematic uncer- 
tainties. All detector components are reproduced in de- 
tail with the GEANT4.9.3.p02 particle transport simu- 
lation code [28, 29]. The simulations additionally uti- 
lize a number of specialized physics packages that are 
tailored to low-energy electromagnetic processes, such as 
G4LowEnergyCompton [30], which take into account mod- 
ifications to the Compton process by the binding and 
kinetic energies of the electrons. 

Since the detector components subtend angles that are 
a non-negligible fraction of the central scattering angles, 
no data set features a perfectly monocncrgctic energy 
deposition, but instead sees a range of energies. The dis- 
tribution of these energies must be understood if the un- 
derlying energy-dependent scintillation light yield is to 
be determined. Figure 4 shows the distribution of raw 
energy deposition for all scattering angles. That is, no 
convolution or energy threshold has been applied, and 
these histograms represent the true distribution of recoil 
energies. Two important features of these distributions 
are apparent. First, the distributions are not symmet- 
ric, and instead feature an extended high-energy "tail". 
Second, the peak of the distributions arc shifted rela- 
tive to the energy values predicted by Eq. (I). These 
effects are a result of the fact that, although the dis- 
tribution of scattering angles is nearly symmetric, the 
recoil energy becomes quadratic in 9 for small 9 (that is, 
E eT « E 2 9 2 /2m e c 2 ). Uniform intervals of energy, AE CI , 
therefore span intervals of 9 that scale as and hence 
smaller angles contribute more to the energy distribution 
than larger angles. We take the central energy to be the 
peak position of the raw MC E CI distribution, and the 
I a range given by a region covering 68.3% of the total 
spectrum, shown in Table I. These bounds are chosen 
in a way such that the differential rate is equal at both 
boundaries (i.e. the dashed-blue lines in Figure 4 cross 



the black histogram at the same vertical position). 

Events in which the 7 ray scatters multiple times in the 
LXe, or in detector materials other than active LXe or 
Nal, represent background populations that can be well 
understood with the MC. The contributions from these 
backgrounds arc small, since the small size of the LXe 
cell makes the former ( "multiple scatters" ) unlikely, and 
the latter ("materials scatters") are reduced by consid- 
ering only events within the full absorption peak in the 
Nal. The simulations indicate that multiple scatters con- 
tribute on average 1.6% to the total signal, and materials 
scatters on average 5.8%. 

The results of the MC simulations are compared to 
data via the likelihood function (see Section IV) in order 
to extract the light-yield values. In an effort to account 
for potential misalignment of experimental components, 
central scattering angles are simulated additionally with 
9 ± 0.125°, and the range of resulting light yield values 
is taken as the geometrical systematic uncertainty. 



IV. ANALYSIS AND RESULTS 

A. Data selection 

As mentioned in Section II, four signals are recorded 
for each trigger: the PMT coupled to the Nal crystal, 
the two PMTs immersed in the LXe, and the output 
of the TAC. Data quality cuts are performed to remove 
events that show significant shift or fluctuations in the 
baseline of the digitized trace. The main data selection 
is performed based upon the Nal and TAC signals. The 
two-dimensional distribution of these parameters, for the 
central scattering angle of 8.5°, is shown in Figure 5. The 
peak along the horizontal axis at roughly 4 ns corresponds 
to the Time of Flight (ToF) of the 7 ray between the LXe 
and Nal. The Compton continuum tends towards higher 
TAC values for lower energy depositions; this 'walk ef- 
fect' is a result of the fact that the start and stop signals 
controlling the TAC are generated by a leading-edge dis- 
criminator with a threshold that is given by an absolute 
voltage amplitude. Therefore, a small signal will trigger 
the discriminator at a later time than a larger — but oth- 
erwise identically shaped — signal. Though this feature 
could be removed by the use of a constant-fraction dis- 
criminator, it is inconsequential to the analysis because 
only events falling within the range [— 1 <7,+3 a] of the 
137 Cs full- absorption peak, located at (662kcV-i? or ), are 
considered. This asymmetric window is chosen in order 
to minimize the contribution from 7 rays that may un- 
dergo a small-angle Compton scatter along the transit 
from the LXe to the Nal. A similar effect occurs with 
LXe trigger signals, but is significantly less pronounced 
than the Nal due to LXe's much faster scintillation pulse 
rise time. A more significant effect on the ToF measure- 
ment for small LXe scintillation signals occurs due to the 
O(10 ns) scintillation emission timescale [32]; for scintil- 
lation signals smaller than ~10 PE the arrival of the first 
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FIG. 4: Monte-Carlo distributions of raw energy deposition of Compton electrons (i.e. no convolution or energy thresholds 
applied) for all central scattering angles. Solid-blue lines indicate the peak position taken as the central energy, and the the 
blue-dashed lines enclose 68.3% of the distribution, taken as the 1 a spread in the energy deposition. The feature in the spectra 
of the top-two plots at ~5keV is due to Xe's L-shells whose binding energies are at 5.45 keV, 5.10 keV, and 4.78 keV [31]. 
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FIG. 5: (color online) Energy deposited in the Nal scintillator 
versus the time of flight of the 7 rays for the data with a 
central scattering angle of 8.5° (E cr = 7.8I4 4). The green 
(solid) box indicates the selection cuts used for the signal, 
blue (dashed) boxes are the sidebands used to estimate the 
background from accidental coincidences. 



PE can be delayed, resulting in a corresponding tail of 
the ToF distribution. For this purpose, asymmetric ToF 
bounds are chosen in order to collect all emitted LXe 
photons. 

The green box in Figure 5 indicates the bounds of 
the event selection cut. Identical peak-selection cuts are 
applied to the MC events. Though the distribution of 
events in ToF is highly peaked, there exists a population 
of events that is uniform in ToF, resulting from acciden- 
tal coincidences between LXe and Nal signals. The con- 
tribution of this background is estimated from the ToF 
sidebands enclosed by the dashed-blue boxes in Figure 5. 
The spectrum of LXe scintillation of the events in these 
sidebands shows no dependence on the ToF, and there- 
fore the sideband boundaries are chosen to maximize the 
statistics of the estimate. 

A potential additional background that is coincident 
with the Nal signal arises from an observed effect 
whereby one of the PMTs gives a signal in response to 
a particle scatter in the PMT alone, likely the dynode 
chain. This effect has been confirmed by placing one of 
the PMTs used here in a black box without any scin- 
tillator and comparing the dark-current spectrum with 
and without a 137 Cs source in the vicinity. A PMT sig- 
nal of this sort is unconnected with any light emission, 
and therefore does not produce signals coincident in the 
two LXe PMTs. In order to reduce this background, 
the 137 Cs collimator is placed as close as possible to the 
LXe cryostat, resulting in a small spot size. Since this 
background is absent from signals that are coincident in 
both LXe PMTs, the data from each scattering angle are 
each separately fit with the hardware- imposed N>1 co- 
incidence requirement, and additionally while imposing 



an N>2 requirement in software. The spectra obtained 
when requiring this more stringent selection requirement 
are corrected with a simulated efficiency curve, and the 
discrepancy between the fit results for N>1 and N>2 
conditions is treated as a systematic uncertainty, shown 
in Table I. 



B. Light-yield determination and results 

For each measured scattering angle, the light yield is 
determined by iteratively transforming the correspond- 
ing MC results into an expected scintillation distribution 
until the likelihood function is maximized. The deposited 
energy, Ei, of each MC event i, is scaled by an energy- 
dependent light yield, LY(Ei), to an expected scintilla- 
tion signal, 



(Si) =EiX LY(Ei), 



(2) 



in units of PE. Because the raw energy distributions are 
not monoenergetic, LY(Ei) should allow for a nonzero 
slope in the region of the peak energy (see Figure 4). 
Outside the peak energy, LY(Ei) should flatten, so that 
the tails of the distributions do not become scaled by 
unphysically large (or small) light yields. Any generic 
sigmoid function can accomplish this purpose; we choose 
here the error function, such that LY(Ei) ~ erf[ln(.Ej)]. 
The argument of the error function is chosen to be log- 
arithmic because the light yield is expected to change 
more at small energies than at large energies, and ad- 
ditionally the raw energy distributions are skewed with 
high-energy tails; a logarithmic energy scale is the sim- 
plest way to capture these features. The LY(Ei) function 
is shifted horizontally so that it is centered at the peak of 
the energy distribution, E c , given a vertical offset, LY3, 
and allowed to vary roughly within the range [E-,E + \ 
(the la bounds of the energy distribution). With these 
adjustments, we take the functional form of the light yield 
to be 



LY[Ei) = LY 



m Er- 



in 



E c 



erf 



\n(E + /E c ) 



(3) 

where E c and E+ are fixed constants for each scattering 
angle, and m and LYq are free parameters. The prefac- 
tors and argument factors are chosen so that this generic 
function has the convenient property that when E\ = E c , 
both the function and its first derivative take on the sim- 
ple forms LY(E C ) = LY and LY'(E C ) = m. The la 
lower bound on the energy distribution, £7_, does not 
explicitly appear in Eq. (3), but is implicitly respected 
because the energy distributions below 10° in Figure 4 
are nearly symmetric in log space. 

The probability to obtain an integer number of PE, j, 
for a given MC event is found by assuming that the (Si) 
from Eq. (2) each define the mean of a Poisson distribu- 
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tion, 



(Si 



-(Si) 



(4) 



where Sij is the probability to obtain j photoclcctrons 
from the i th MC event. The expected number of events 
with j photoelectrons for the entire MC run, Sj, is 
then given by summing the Sij over all MC events, 
Sj ~J2i Sij ■ 

The spectrum Sj is convolved with the measured 
single-PE response of the LXe PMTs to produce an ex- 
pected measured signal (note that a measured signal can 
have fractional number of PE). The top and bottom 
PMTs have single-PE resolutions (<r//i) of 54% an 55%, 
respectively. An additional Gaussian convolution is ap- 
plied to account for extra fluctuations that might arise 
from, for example, effects such as electron-ion recombi- 
nation or position-dependent geometrical light collection. 
In order to avoid making assumptions about the behavior 
of these effects, the functional form of the Gaussian con- 
volution's width, to, is taken to be a power law in j with 
a free exponent, w(j) = w\j W2 + u>3, which adds three 
free nuisance parameters to the fit (u>i, u>2, W3). Finally, 
this spectrum is re-binned to match the bin size of the 
real data, and then normalized according to the source 
activity, A. This normalization is treated as a free param- 
eter constrained with a Gaussian prior to account for the 
uncertainty in the activity. This uncertainty is found by 
leaving A completely free in initial fits to the 8.5°, 16.25°, 
and 34.5° datasets, and combining the resulting statisti- 
cal uncertainties to form A's prior. The spectrum of acci- 
dentals background (estimated from the ToF sidebands) 
is added to the MC spectrum, forming the expected con- 
tent of the k th bin, , which is implicitly a function of the 
six fit parameters, vu = i/],(LYo,m,Wi, 11)2,11)3, A). The 
standard Poisson log-likelihood function, weighted with 
a Gaussian prior for the source activity, can be written 
as 



= 7~A n fc mi/ fc 



Vk) 



A 2 - 2Afi A 



26/ 



(5) 



where is the histogram of real data, [ia and 8a are 
the mean and width, respectively, of the Gaussian prior 
constraining the source activity. The sum is carried out 
over bins within the fit range (f.r.), which is chosen to 
roughly cover the peak of the observed spectrum (be- 
cause LY(E) is featureless outside this peak region) and 
to avoid contributions from the long tails of the energy 
distributions. In order to understand how the choice of 
the f.r. affects the resulting LYo, each f.r. is later varied 
and the observed discrepancies in LYo are included as a 
systematic uncertainty (Section IV C). 

The likelihood function (Eq. (5)) is proportional to the 
Bayesian posterior probability density function (PDF). 
This PDF is sampled with the Metropolis-Hastings algo- 
rithm [33, 34], which is an iterative Markov-Chain pro- 
cess that is well suited for producing a random sample of 



data from a multidimensional PDF for which evaluation 
may be computationally intensive. It has the additional 
useful property that, because the output is a set of data, 
marginalizing the posterior PDF over nuisance parame- 
ters (m, Wi, and A) is as simple as histogramming the de- 
sired parameter, LYq. The data collected at all scattering 
angles arc shown in Figure 6, along with the correspond- 
ing best-fit MC spectra. The statistical uncertainties are 
given by the la contour of the marginal posterior PDF 
in LYq. The subtle L-shell feature seen in the raw-energy 
spectra (Figure 4) at 4.25° and 5.25° is not visible in the 
scintillation data due to Poisson smearing. 

Following the procedure of [14], we present lZ e , the LY 
results relative to the 32.1 keV emission of 83m Kr, shown 
in Table I and Figure 7 (the reason for this normalization 
is discussed in Section VB). Also shown are the corre- 
sponding Ti e values from the 57 Co and 83m Kr sources, 
along with the results of Obodovskii and Ospanov us- 
ing induced X-rays [35], the recent measurement by 
Aprile et al. [14], and the prediction from the Noble Ele- 
ment Simulation Technique (NEST) version 0.98 [36, 37], 
which is a constrained implementation of the Thomas- 
Imel electron- ion recombination model [38]. The values 
from Obodovskii and the curve from NEST have been 
vertically normalized such that their interpolated values 
at 32.1 keV arc unity. The gray band indicates the la 
allowed LY models used in energy-threshold determina- 
tions, discussed in Section VC. 



C. Systematic uncertainties 

Five systematic uncertainties are considered, as shown 
in Table I. 

• <jW results from potential misalignment of the ex- 
perimental components (discussed in Section III), 
and is studied by simulating central angles with a 
shift of ±0.125° and determining the best LYo as 
before. 

• a^ 2 ' , the uncertainty from the choice of fit ranges, is 
studied by modulating the chosen fit ranges at the 
level of 20% and taking the observed differences in 
the best-fit LYq values. The spectra from 6.25° and 
8.5° in Figure 6 show a discrepancy between data 
and MC in the high-energy tail; this discrepancy 
is absent in the fit obtained with the alternate fit 
range, with little effect on the reconstructed light 
yield. 

• a^ indicates the dependence of the best-fit LYo on 
the source-activity parameter, A. This systematic 
is determined by using the covariance between LYq 
and A from the fit to estimate how much LYq can 
vary within the allowed uncertainty on A, given by 



(3) COY (LYq, A) 
a y ' = 5 6a, 
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FIG. 6: (color online) Comparison between the scintillation spectra of real data (red boxes) and the simulated data using the 
best fit parameters (black) for all scattering angles. Blue-dashed lines indicate the main fit range. The discrepancy observed 
in the high-energy tails of the data obtained at 6.25° and 8.5° is absent when alternate fit ranges are used, with little effect on 
the reconstructed light yield (see text). 
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TABLE I: Results of the light-yield measurements. 9 C is the central angle of the dataset; E CI is the central energy of the energy 
distribution; TZ e is the zero-field central relative light yield value (relative to the scintillation emission at 32.1 keV); <r st is the 
statistical uncertainty; o~s y l is the systematic uncertainty resulting from potential misalignment of experimental components; 

(2) (3) 

CTsys is the systematic uncertainty associated with the choice of fit range; o-gys is the systematic uncertainty associated with source 
activity; a^yl indicates the discrepancy introduced between 1-fold and 2-fold coincidence requirements on the LXe PMTs; an 
additional systematic uncertainty of 1.5% is applicable to all values in the third column, which arises from variations in results 
of weekly 57 Co calibrations. g(450) is the scintillation quenching factor at an applied field of 450 V/cm; the first uncertainties 
are statistical, the second systematic. 
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FIG. 7: (color online) Results of the light yield relative to 
that of the 32.1 keV emission of 83m Kr, TZ e . The current 
work (red) shows statistical uncertainties as vertical lines, sys- 
tematic uncertainties as light, shaded rectangles, and the 1 a 
spread in the distribution of electron recoil energies as hori- 
zontal lines. Also shown are the results from studies with X- 
rays [35] (blue), the recent Compton-scatter study by Aprile 
et al. [14] (purple) and the model prediction of NEST [36, 37] 
(green). The gray band indicates the la range of 1Z £ models 
used to determine the energy thresholds of four recent LXe 
dark-matter searches. 



where 8a is the uncertainty in the source activity 
(as in Eq. (5)) and u\ is the variance of A from the 
fit. The factor cov{LYq, A) / a\ gives the slope of 
LYq versus A. 

• a' 4 ' quantifies the uncertainty associated with the 
choice of the PMT coincidence requirement. An 
N = 2 coincidence requirement on the two LXe 
PMTs is separately imposed, correcting the result- 
ing scintillation spectrum by a simulated coinci- 
dence efficiency curve, and performing the fits again 
for LY . 

• cr^ 5 ) is a 1.5% relative systematic from fluctuations 
in the PMT gains and weekly 57 Co calibrations. 

These systematic uncertainties are combined in quadra- 
ture to form the systematic error bars in Figure 7, and 
the first four are shown in Table I. In the lowest energy, 
the dominating systematic is er' 4 ' with a contribution of 
38%; this systematic rapidly decreases to 1% by 8.5° and 
zero beyond. 

D. Field dependence 

The previous results all pertain to the light yield of 
LXe with no applied electric fields. As mentioned in Sec- 
tion II, data were also collected with an applied field 
of 450 V/cm for a subset of scattering angles in order to 
study the scintillation quenching of LXe at the lowest en- 
ergies. The data collected with this field are fit using the 
same procedure as before, resulting in a set of posterior 
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FIG. 8: (color online) The quenching of the scintillation signal 
with an applied electric field of 450V/cm. Vertical lines rep- 
resent statistical uncertainties, grey bars represent systematic 
uncertainties, and horizontal lines are the 1 a spread in the 
distribution of electron recoil energies. Also show are the pa- 
rameterized predictions from [13] (blue circles) and 57 Co field 
quenching [39] (purple diamonds) at 400 V/cm and 500 V/cm. 
The prediction of the NEST model [36, 37] for quenching at 
450 V/cm is indicated by the green curve. 



PDFs for the light yield. The last row of Figure 6 shows 
the measured and best-fit spectra of the three scattering 
angles collected. These PDFs arc convolved with their 
corresponding zero-field light yield PDFs to obtain pos- 
terior PDFs of their ratio, known as the field-quenching 
value, g(450), shown in Table I. For each scattering an- 
gle with applied field, the 450 V/cm data and the zero- 
field data were taken consecutively. Therefore, any po- 
tential misalignment of experimental components will be 
unrelated to the applied field. The resulting scintilla- 
tion quenching values, along with those simultaneously 
obtained for 57 Co and 83m Kr, are shown in Figure 8. 
Also shown is the predicted scintillation quenching of the 
NEST model. 



support an average value of g(450) = 0.74±0.11. For the 
NEST prediction of this quantity shown in Figure 8, the 
horizontal scale indicates the energy of the primary 7 ray 
(not electronic-recoil energy), and is therefore in princi- 
ple distinct from Compton scatters. The feature in the 
NEST curve between ~15 keV and ^50 keV is an indirect 
result of photoabsorption on if -shell electrons, and would 
be absent for Compton scatters of this energy. However, 
the distinction between Compton scatters and photoab- 
sorptions disappears at low energies [36, 40], where the 
recombination probability becomes independent of stop- 
ping power, and instead depends only on the total num- 
ber of charges produced. It is therefore an applicable 
prediction of our results in this energy regime. 

It is interesting to note that the data obtained from 
X-rays [35] show an increased light yield at 7.84 keV 
compared with the data obtained here from Compton 
scatters, when normalizing their interpolated value at 
32.1 keV. The photoabsorption process that the X-rays 
undergo favors inner-shell electrons (when accessible) 
[41], which means that the recoiling electrons can have 
significantly less energy than the incoming photons be- 
cause they must overcome large binding energies. On 
the other hand, Compton scattering on inner-shell elec- 
trons is suppressed for scattering angles below ^60° [42] . 
Therefore, the two results actually probe LXc's response 
at slightly different electron energies. In principle, the 
axioelectric effect, which has been induced as a possi- 
ble explanation of the observed DAMA annual modula- 
tion signal, would be similar to the photoelectric effect. 
However there is of course an overlap of effects, since 
low-energy Compton scatters do also probe inner-shell 
electrons, as can be seen by the L-shell feature in Figure 
4. 

The data reported by Aprile et al. [14] show good 
agreement with the present results above ^10keV, but 
show a separation below this energy. Considering both 
statistical and systematic uncertainties gives a maximum 
discrepancy of 1.7cr at ^5keV and 1.4cr at ~1.5keV. 



B. The 9.4 keV anomaly 



V. DISCUSSION 

A. Comparison of results 

The results presented here represent the first obser- 
vation of LXe scintillation light from electronic recoils 
down to 1.5 keV, and additionally measure the behav- 
ior of this scintillation emission under the application of 
a static electric field. The general behavior — that of re- 
duced LY for decreasing energies — is predicted by a num- 
ber of methods (see [36] and references therein), and is 
understood as being due to reduced electron-ion recom- 
bination. Below lOkeV, the data show no significant en- 
ergy dependence on the strength of field quenching, but 



The discrepancy seen in the LY of the 9.4 keV emission 
from 83m Kr deserves attention. The energy of this decay 
is carried mostly by internal conversion electrons emitted 
from the inner shell [43], however, this data point is in- 
consistent also with the X-ray data, for which the process 
should in principle be similar. One notable characteris- 
tic of the 9.4 keV emission is that it quickly follows the 
32.1 keV emission of the same nucleus, with a half-life of 
154.4ns [44]. It was pointed out by [45] that the 32.1kcV 
emission could leave behind a cloud of electron-ion pairs, 
close to the mother nucleus, that fail to recombine. The 
electrons (ions) produced by the 9.4 keV emission could 
then potentially have an additional supply of left-over 
ions (electrons) with which to recombine, producing more 
scintillation photons than would be observed normally. 
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We test this possibility by investigating the time de- 
pendence of the 9.4 keV LY, since the scintillation en- 
hancement should disappear as these charges diffuse away 
over time. Figure 9 (top) shows the size of the 9.4 keV 
scintillation signal at zero field as a function of delay be- 
tween 32.1 keV and 9.4 keV emissions, At. A clear rise 
in scintillation signal is seen for delay times <400ns. Al- 
though making a precise prediction for the diffusion rate 
of the electrons and ions left over from the 32.1 keV emis- 
sion is difficult, a simple order-of-magnitude estimate can 
provide a useful expectation. For an initial population of 
particles distributed over a volume of space, a charac- 
teristic timescale for the diffusion process, td, can be 
estimated through dimensional analysis by 



TD= D> 



(7) 



where a is a length characterizing the size of the initial 
particle distribution, and D is the diffusion coefficient of 
those particles. The typical range of a 30keV electron 
is on the order of ^10 fim in LXe [46] and the diffusion 
coefficient of electrons is D_ ~ 60cm 2 /s [47]. The ionic 
mobility of Xe + is given in the literature as /i+ ps 4 x 
10~ 3 cm 2 /Vs [48], which can be connected to its diffusion 
coefficient by the Nernst-Einstein relation [49], 



D+ = — M+, 

e 



(8) 



giving D + ~6x lCT 5 cm 2 /s at T = 180 K. With these 
values for a and D±, we can estimate the diffusion 
timescales as 



T D _ = O(10 ns) 
t d+ = O(10 ms) 



(e") 
(Xe + ) 



(9) 
(10) 



If the initial particle distribution is taken to be Gaussian, 
it is straightforward to show that the central number den- 
sity, n c , evolves with time according to 



oc (2t + r D ) 



-3/2 



(11) 



To illustrate how this compares to the observed behav- 
ior, Eq. (11) is plotted in Figure 9 (top) for m = 10 ns 
(dashed-cyan curve). The same curve using td = 10ms 
appears simply as a straight line on this scale, which is 
not surprising — given the 154 ns half-life of the 9.4 keV 
state, it would be unlikely to see a change in its LY due 
to Xe + diffusion. However, these leftover ions could still 
enhance the scintillation signal of the 9.4 keV transition, 
albeit with no observable time dependence. The electron 
diffusion timescale estimated here is consistent with the 
observed time dependence. It is therefore reasonable to 
conclude that the anomalous LY of the 9.4 keV emission 
from 83m Kr is due to left-over charges from the preceding 
32.1 keV emission. 

An applied electric field provides an extra means (in 
addition to diffusion) by which electrons can leave the 
vicinity of the decaying nucleus. The timescale at which 
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FIG. 9: (color online) The distribution of scintillation sig- 
nals from the zero-field 9.4 keV emission of 83m Kr (top) as a 
function of the time since the same nucleus's 32.1 keV transi- 
tion, At. The dashed line indicates the expected time depen- 
dence for diffusion of electrons from a 32.1 keV interaction. 
Also shown are the corresponding distributions for 9.4 keV at 
450 V/cm (middle) and 32.1 keV at zero field (bottom). In all 
panels, the blue histograms indicate the mean of the scintil- 
lation distributions in progressive time slices. 



this occurs can be estimated by the time required for the 
charges to traverse the affected region, given by 



— = — = 5 ns, 

|E|/x_ v d 



(12) 



where /i_ is the electron mobility, |E| is the electric field 
strength, and Vd « 2mm//is is the drift velocity of elec- 
trons [50] . Given that this process is of similar order to 
the electron diffusion process, one expects the scintilla- 
tion signal at |E| = 450 V/cm to show less of an effect 
from electron diffusion. This time dependence is shown 
in Figure 9 (middle), and exhibits a reduced scintillation 
increase at low At. 

Though the average zero-field LY of these data gives 
lZ e = 1.10 as quoted in Table I, this value depends 
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on the range of At considered. Since the characteris- 
tic PMT output from a 83m Kr decay is a trace contain- 
ing two scintillation pulses, the pulse-finding algorithm 
used must correctly identify this feature. However, this 
can be difficult for small At, when the pulses begin to 
overlap. The efficiency for the pulse-finding algorithm 
implemented here to correctly identify such double-pulse 
events is unity for At > 150 ns, but this would not 
necessarily be true for other detectors utilizing a dif- 
ferent analog bandwidth, data acquisition system, and 
data processing techniques. For these reasons, 83m Kr's 
9.4 keV transition is not very well suited as a "standard 
candle" calibration source in the way that 57 Co is often 
implemented. In principle, the 9.4 keV transition could 
be used for standardized calibrations if only events in 
which At > 400 ns are considered. However, although 
the results presented here show no significant time de- 
pendence above this value, the results in [14] do exhibit 
a continued decrease until roughly 1 /xs, and therefore one 
cannot predict how significant this observed time depen- 
dence might be in various other systems. On the other 
hand, the 32.1 keV transition shows no dependence on 
At, as expected, which we show here in Figure 9 (bot- 
tom). We therefore conclude, in agreement with [14], 
that the 32.1 keV transition of 83m Kr provides a good 
calibration source by which to compare the scintillation 
response of various LXe detectors. 



C. Impact on dark matter searches 

The primary motivation for the present study is to 
learn whether existing LXe dark matter search results 
have energy thresholds that are low enough to probe a 
2-5 keV electronic-recoil peak, as potentially observed in 
the DAMA/LIBRA experiment. To illustrate how our 
results can address this question, we consider four ex- 
isting LXe dark-matter-search results: ZEPLIN-III [51], 
XENON10 [3], XENON100 [52], and XMASS [7]. We 
wish to determine their electronic-recoil energy thresh- 
olds based on their quoted scintillation thresholds. Given 
an electronic recoil that deposits energy E cv , the average 
scintillation signal, 51, in units of PE will be given by 



51 = E CI x f Co (E er ) x LY Co x 



g(|E|) 

qco 



(13) 



where /co(-Ecr) is the ratio of the zero-field LY at E er 
to that from 57 Co, E is the applied field, LYqo is the 
57 Co light yield (in PE/kcV) at E, qc a is the scintillation 
quenching of 57 Co at E, and g(|E|) is the scintillation 
quenching at energy E er and field E. To determine an 
energy threshold, Eq. (13) is inverted and evaluated for 
the quoted scintillation threshold, 51thr- 

All four experiments utilize different applied elec- 
tric fields, and therefore we must extrapolate the field- 
quenching results reported here to the appropriate val- 
ues. This extrapolation introduces an uncertainty in the 



calculated energy threshold that ranges from nonexis- 
tent in XMASS (|E| = 0) to negligible in XENON100 
(|E| = 530V/cm) to considerable in ZEPLIN-III (|E| = 
3400 V/cm). In order to do so, we use an empiri- 
cal parameterization of the field quenching inspired by 
the Thomas-Imel electron-ion recombination model [38] , 
which was applied to LXe's 83m Kr response in [13] as 



g(|E|) = oia 3 |E| In 1 



1 



o 3 |E| 



(14) 



where a\ and ai are free parameters, with a\ describing 
the overall strength of the field quenching, and ai de- 
scribing the field dependence of this quenching. In [13], 
it was found that the energy dependence of 02 is much 
less significant than that of a±. Therefore, we use here 
a-2 = (8.3 ± 1.7) x 10~ 4 cm/V, which is the average of 
the values for 9.4keV and 32.1keV reported in [13]. The 
value of ai is chosen so that the <z(|E|) function is consis- 
tent with our g(450) value, which is taken from the aver- 
age of data collected at 4.25°, 5.25°, and 8.5° (indicated 
by the red circle in Figure 10). Combining the uncertain- 
ties in ai and g(450) produces the bands show in Figure 
10, which is taken to represent the energy-averaged field 
quenching below ~T0kcV. 
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FIG. 10: (color online) The 1<t and 2a bands of the scin- 
tillation field quenching below ~10keV, g(|E|), used in the 
determination of -Ethr. Also indicated are the fields utilized 
by the four dark-matter experiments considered in the text. 
The red circle indicates the measured g(450), averaged from 
data obtained at 4.25°, 5.25°, and 8.5°. 



Also needed in the determination of the electronic- 
recoil energy threshold, Ethr, is a model of /co(-Eer). For 
this, a range of models are taken that fit our lZ e data, and 
whose la span is indicated by the gray band in Figure 7. 
The uncertainties on these three parameters (a 2 , <;(450), 
and fco(E cr )) are convolved to produce likelihood curves 
for the resulting E t hr of the four experiments considered 
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here. The results are shown in Table II. It is clear that 
all four experiments, even in the presence of the sharply 
falling 7Z e observed here, have sensitivity to all or part 
of the 2-5 keV range favored by the DAMA results. 

TABLE II: Four recent dark-matter searches using LXe: the 
second science run of ZEPLIN-III [51], results of XENON10 
[3], the recent 225 live days reported from XENON100 [52], 
and the results of XMASS [7] . Shown are the applied electric 
fields used by each (|E|), their quoted scintillation thresholds 
(Slthr), their 57 Co light yield (LYco), and their electronic- 
recoil energy thresholds using this work (Ethr)- 



Experiment 


|E| (V/cm) 


Slthr (PE) 




E thT (keV) 


ZEPLIN-III 


3400 


2.6 


1.3 


2 4 +o.5 


XENON10 


730 


4.4 


3.0 




XENON100 


530 


3.0 


2.3 


1 7 +0A 

L -> -0.3 


XMASS 





4.0 


14.7 


1 1+0.4 



VI. SUMMARY 



roughly 40% of its value at higher energies. With the ap- 
plication of a static electric field of 450 V/cm, we observe 
a reduction of the scintillation signal of roughly 75% rel- 
ative to the value at zero field, and see no significant 
energy dependence on this value between 1.5 keV and 
7.8 keV. With these values, we are able to extrapolate 
the electronic-recoil energy thresholds of the ZEPLIN- 
III [51], XENON10 [3], XENON100 [52], and XMASS 
[7] experiments. These experiments report scintillation 
thresholds of 2.6 PE, 4.4 PE, 3.0 PE, and 4.0 PE, which, 
when applied with the results presented here, give energy 
thresholds of 2.4keV, 1.8 keV, 1.7 keV, and 1.1 keV, re- 
spectively. We additionally investigate a discrepancy be- 
tween the LY from the 9.4 keV emission of 83m Kr (which 
has in the past been considered for use as a standard 
calibration source) and other observed LY values nearby 
in energy. We observe a time dependence of this scintil- 
lation signal which we show is most likely due to diffu- 
sion of electrons from the preceding 32.1 keV decay of the 
same nucleus. From this observation, we conclude that 
the 9.4 keV peak is not an optimum standard calibration 
feature, and instead advocate the use of the 32.1 keV for 
this purpose. 

Acknowledgments 



The work presented here details a study of LXe's scin- 
tillation response to electronic recoils as low as 1.5 keV. 
The proportionality between deposited energy and scin- 
tillation signal, or "light yield" (LY), is observed to drop 
with decreasing energy beginning at ~10keV, to a level 



This work was supported by the Swiss National Foun- 
dation grant No. 200020-138225. A.K. and A.M. received 
support from the University of Zurich Forschungskredit 
fellowship, Nos. 57161702 and 57161703, respectively. 



[1] G. Bertone, D. Hooper, and J. Silk, Phys. Rept. 405, 279 

(2005), hep-ph/0404175. 
[2] E. Aprile and T. Doke, Rev. Mod. Phys. 82, 2053 (2010), 

0910.4956. 

[3] J. Angle et al. (XENON10 Collaboration), Phys. Rev. 

Lett. 100, 021303 (2008), 0706.0039. 
[4] V. Lebedenko, H. Araujo, E. Barnes, A. Bewick, R. Cash- 
more, et al., Phys. Rev. D80, 052010 (2009), 0812.1150. 
[5] E. Aprile et al. (XENON100 Collaboration), Phys. Rev. 

Lett. 107, 131302 (2011), 1104.2549. 
[6] D. Akerib, X. Bai, S. Bedikian, E. Bernard, A. Bernstein, 

et al., Nucl. Instrum. Meth. A668, 1 (2012), 1108.1836. 
[7] K. Abe et al. (XMASS Collaboration), Phys.Lett. B719, 

78 (2013), 1211.5404. 
[8] V. Chepel et al., Astropart. Phys. 26, 58 (2006). 
[9] E. Aprile, L. Baudis, B. Choi, K. Giboni, K. Lim, et al., 

Phys. Rev. C79, 045807 (2009), 0810.0274. 
[10] A. Manzur, A. Curioni, L. Kastens, D. McKinsey, K. Ni, 

et al., Phys. Rev. C81, 025808 (2010), 0909.1063. 
[11] P. Sorensen and C. E. Dahl, Phys. Rev. D83, 063501 

(2011), 1101.6080. 
[12] G. Plante, E. Aprile, R. Budnik, B. Choi, K. Giboni, 

et al., Phys. Rev. C84, 045805 (2011), 1104.2587. 
[13] A. Manalaysay, T. Marrodan Undagoitia, A. Askin, 

L. Baudis, A. Behrens, et al., Rev. Sci. Instrum. 81, 

073303 (2010), 0908.0616. 



[14] E. Aprile, R. Budnik, B. Choi, H. Contreras, K. Giboni, 
et al., Phys. Rev. D86, 112004 (2012), 1209.3658. 

[15] F. T. Avignone III, R. L. Brodzinski, S. Dimopou- 
los, G. D. Starkman, A. K. Drukier, D. N. Spergel, 
G. Gelmini, and B. W. Lynn, Phys. Rev. D 35, 
2752 (1987), URL http://link.aps.org/doi/10.1103/ 
PhysRevD.35.2752. 

[16] C. Aalseth et al. (CoGeNT collaboration), Phys. Rev. 
Lett. 106, 131301 (2011), 1002.4703. 

[17] R. Bernabei et al. (DAMA Collaboration, LIBRA Col- 
laboration), Eur. Phys. J. C67, 39 (2010), 1002.1028. 

[18] R. Bernabei, P. Belli, F. Montecchia, F. Nozzoli, F. Cap- 
pella, et al., Int. J. Mod. Phys. A21, 1445 (2006), astro- 
ph/0511262. 

[19] Z. Ahmed et al. (CDMS Collaboration), Phys. Rev. Lett. 

103, 141802 (2009), 0902.4693. 
[20] A. Manalaysay, Ph.D. thesis, University of Florida, 

Gainesville, Florida (2009). 
[21] M. Yamashita et al., Nucl. Instrum. Meth. A 535, 692 

(2004), ISSN 0168-9002. 
[22] A. I. Bolozdynya, Nucl. Instrum. and Meth. A 422, 314 

(1999). 

[23] R. Dossi, A. Ianni, G. Ranucci, and O. Y. Smirnov, Nucl. 

Instrum. Meth. A451, 623 (2000). 
[24] E. Aprile et al. (XENON100 Collaboration), Astropart. 

Phys. 35, 573 (2012), 1107.2155. 



14 



[25] L. Kastens, S. Cahn, A. Manzur, and D. McKinsey, Phys. 
Rev. C80, 045809 (2009), 0905.1766. 

[26] COMSOL Multiphysics® , http : //www. comsol . com. 

[27] H. Dujmovic, Simulation and optimization of the elec- 
tric field in a liquid xenon time projection chamber 
(2012), Bachelor's Thesis, University of Zurich, URL 
http : //www.physik .uzh. ch/groups/groupbaudis/ 
darkmatter/ theses/ xenon/bachelor_du jmovic . pdf . 

[28] S. Agostinelli et al. (GEANT4 Collaboration), Nucl. In- 
strum. Meth. A506, 250 (2003). 

[29] J. Allison, K. Amako, J. Apostolakis, H. Araujo, 
P. Dubois, et al., IEEE Trans. Nucl. Sci. 53, 270 (2006). 

[30] F. Longo, L. Pandola, and M. Pia, in Nuclear Science 
Symposium Conference Record, 2008. NSS '08. IEEE 
(Oct.), pp. 2865-2868, ISSN 1095-7863. 

[31] R. Deslattes, E. Kessler, P. Indelicato, L. de Billy, E. Lin- 
droth, J. Anton, J. Coursey, D. Schwab, C. Chang, 
R. Sukumar, et al., X-ray Transition Energies (ver- 
sion 1.2) (National Institutes of Standards and Tech- 
nology, Gaithersburg, MD, 2005), [Online] Available: 
http://physics.nist.gov/XrayTrans [2013, March 11]. 

[32] A. Hitachi, T. Takahashi, N. Funayama, K. Ma- 
suda, J. Kikuchi, and T. Doke, Phys. Rev. B 27, 
5279 (1983), URL http://link.aps.org/doi/10.1103/ 
PhysRevB.27.5279. 

[33] S. Chib and E. Greenberg, The American Statistician 49, 
327 (1995). 

[34] P. Saha, Principles of Data Analysis (Capella Archive, 

Great Malvern, UK, 2002). 
[35] I. Obodovskii and K. Ospanov, Instrum. Exp. Tech. 37, 

42 (1994). 

[36] M. Szydagis, N. Barry, K. Kazkaz, J. Mock, D. Stolp, 
et al., JINST 6, P10002 (2011), 1106.1613. 



[37] NEST, URL http://nest.physics.ucdavis.edu. 

[38] J. Thomas and D. A. Imel, Phys. Rev. A 36, 614 (1987). 

[39] E. Aprile, C. Dahl, L. DeViveiros, R. Gaitskell, K. Gi- 

boni, et al., Phys. Rev. Lett. 97, 081302 (2006), astro- 

ph/0601552. 

[40] M. Szydagis (2013), private communication. 

[41] G. F. Knoll, Radiation Detection and Measurement 
(John Wiley Sons; New York, USA, 2000), 3rd ed. 

[42] B. Talukdar, S. Mukhopadhyaya, and D. Chattarji, Phys. 
Lett. A 35, 37 (1971). 

[43] A. Chan et al. (DELPHI Collaboration), IEEE Trans. 
Nucl. Sci. 42, 491 (1995). 

[44] S. Chu, L. Ekstrom, and R. Firestone, WWW Table of 
Radioactive Isotopes, http://ie.lbl.gov/toi. 

[45] G. Plante (2012), private communication. 

[46] M. Berger, J. Coursey, M. Zucker, and J. Chang, ESTAR, 
PSTAR, AND ASTAR: Computer Programs for Calcu- 
lating Stopping-Power and Range Tables for Electrons, 
Protons, and Helium Ions (National Institutes of Stan- 
dards and Technology, Gaithersburg, MD, 2005), version 
1.2.3, URL http://physics.nist.gov/Star. 

[47] T. Doke, Nucl. Instrum. Meth. A 196, 87 (1982). 

[48] O. Hilt and W. Schmidt, J. Phys.: Condens. Matter 6, 
L735 (1994). 

[49] O. Hilt, W. Schmidt, and A. Khrapak, IEEE Trans, on 

Dielectrics and Insulation 1, 648 (1994). 
[50] L. S. Miller, S. Howe, and W. E. Spear, Phys. Rev. 166, 

871 (1968). 

[51] D. Y. Akimov, H. Araujo, E. Barnes, V. Belov, A. Be- 
wick, et al., Phys. Lett. B709, 14 (2012), 1110.4769. 

[52] E. Aprile et al. (XENON100 Collaboration), Phys. Rev. 
Lett. 109, 181301 (2012), 1207.5988. 



